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Abstract 

In this paper we study the time evolution of (Renyi) entanglement entropies 
for locally excited states in four dimensional free massless fermionic field theory. 
Locally excited states are defined by being acted by various local operators on 
the ground state. Their excesses are defined by subtracting (Renyi) entanglement 
entropy for the ground state from those for locally excited states. They finally 
approach some constant if the subsystem is given by half of the total space. They 
have spin dependence. They can be interpreted in terms of quasi-particles. 


1 Introduction 


Quantum entanglement is an essential idea that distinguishes quantum physics from clas¬ 
sical physics. There has been a lot of work done to investigate various subfields of physics 
from a perspective of quantum entanglement. For instance, the universal properties of 
quantum entanglement are used to characterize conformal field theories m 121 El S] , as 
well as topologically ordered phases The quantum entanglement also plays an es¬ 

sential role to understand the Hilbert space of gravity mi- While the most commonly 
used quantities for characterizing states are correlation functions of operators, quantum 
entanglement focuses more directly on the Hilbert space. 

Entanglement entropy is defined as the von Neumann entropy of a reduced density 
matrix pa for a subsystem A. It has been observed that the dominant contribution of 
entanglement entropy in ground states of gapped systems comes from near boundary re¬ 
gion of A. This leads to the area law of entanglement entropy [9]. For ground states 
of gapless systems, the contribution of the entanglement from a long distance does not 
necessarily become subdominant. Indeed, in 2d conformal field theories (CFTs), entan¬ 
glement entropy behaves universally as | log l where c is the central charge and l is the size 
of the subsystem A, suggesting that entanglement of any distance equally contributes to 
the entropy. Recently some entanglement measures are introduced to describe the more 
detail structure of the entanglement entropy mm 

The above story is only for static (ground) states and much less is understood for 
the dynamical aspects of the entanglement entropy as well as the entanglement entropy 
for excited states. There are interesting questions about the dynamics of quantum en¬ 
tanglement, such as how the quantum entanglement is created and propagates, how they 
distributes in space time. Those are essential for understanding far-equilibrium states and 
their thermalization [12] , how efficiently quantum dynamics can be simulated in classical 
computers, and where and how the black hole information goes. 

In this paper, we study the time evolution of (Renyi) entanglement entropies by be¬ 
ing acted by a fermionic local operator on the ground state. First we explain quantum 
quenches which are a useful protocol to investigate the dynamical aspects of the entangle¬ 
ment. Suppose a system is prepared for a Hamiltonian H(Xq) which has an experimentally 
controllable parameter A 0 . The state prepared for this system is given by a ground state. 
Then at a certain time the parameter is shifted from Ao to Ai. The prepared state is no 
longer a ground state of H{ Ai) and the state undergoes a time evolution. This kind of 
control is indeed possible in experiments such as cold atom systems m- When parame¬ 
ters are changed globally, these protocols are called as global quenches mt On the other 
hand, it is called as local quenches if parameters are changed locally [1 5] . 

Our protocol is similar to local quenches. In our setup, parameters of Hamiltonian 
are not changed. Instead of being changed them, a local operator O(to,x o) acts on 
the ground state at certain time and creates entangled quasi-particle. In CFTs they 
propagate spherically at the speed of light. If this local operator is inserted outside of 
A, a part of quasi-particles eventually enters into the subsystem A and the rest part 
stays outside. Since they are entangled, it may create an entanglement between the 
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subsystem A and the rest of the system. Indeed this kind of behavior, the increase of 
the entanglement and the saturation, has been observed in various quantum held theories 
[sum mum ism ISM .On the other hand, in holographic held theories the excesses of 
(Renyi) entanglement entropies do not saturate and keep to increase logarithmically even 
if time passes efficiently [191 [231 23J. 

The standard method for computing the entanglement entropy for ground states or 
thermal states in path integral is the replica method: we compute the (Renyi) entangle¬ 
ment entropy S = ^- 0 - log tr^pd first and then take n —* 1 limit, which gives entan¬ 
glement entropy Sa = — ^aPa log pa- Here pa = tr#p and we trace out the degrees of 
freedom outside A (the region B). 

We consider free massless fermionic held theory in 4 dimensional spacetime and take a 
half space as the subsystem A. We generate an excitation by being acted by local operator 
O at a distance l away from the entangling surface and time —t. The state prepared is 
given by a locally excited state. 


|U=VO(-(,-i,x)|0), (1) 

where x = (x 2 ,xs). We define the excesses of (Renyi) entanglement entropies Aby 
subtracting (Renyi) entanglement entropies for the ground state from those for locally 
excited states. 

They do not change for t < l, as expected from the causality. The main difference 
is the hnal values of A S^: In the case of free massless fermionic held case, Ahave 
spin dependence. Their density matrices can depend on the direction of spin because 
the probability with which (anti-)particles are included in A can depend on the direction 
of spin. In the free massless scalar held theory, of course Ado not have such a 
dependence. 

This paper is organized as follows. 

In Sec.2, we explain our setup. In Sec.3, we explain the replica method and the 
analytic continuation which we perform. In Sec.4 we derive the Green function in 4 d 
free massless fermionic held theory on the replica space. In Sec.5 we compute the time 
evolution of A S^\ In Sec .6 we explain the same thing in terms of quasi-particles. In 
Sec.7, we conclude and discuss our results. 

2 Setup 

We study the time evolution of excesses of (Renyi) entanglement entropies for locally 
excited states which are dehned by being acted by various local operators on the ground 
states in following setup. We consider the 4 dimensional free massless fermionic theory, 

Sfermion = ~ J <f 8^, (2) 

where 7 ^ = { 7 *, 7 1 , 7 2 , 7 3 } and t/S = 
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Figure. 1: The location of an local operator and the subsystem A in the Minkowski 
spacetime. 

A local operator O acts on the ground state as in Figure.1. A locally excited state is 
given by 

|#)=jVO(-t,-(,x)|0>, (3) 

where x = ( x 2 ,x 3 ) and A/" is a normalization constant. 

The subsystem A is given by a half of the total space, x 1 > 0 as in Figure. 1. We trace 
out the degrees of freedom in the complement space B outside A and define a reduced 
density matrix, 

pT = tr bP EX (4) 

where p = |T) (T|. 

By using this reduced density matrix, (Renyi) entanglement entropy for locally excited 
state is defined by 

s{ a BX - T^ 1 ^Ap E a X TI ( 5 ) 

By using the reduced density matrix for the ground state = £rg |0) (0|, (Renyi) entan¬ 
glement entropy for the ground state is defined by 

Sa )C = log [tr A ( p f ) n ]. (6) 

We define the excess of (Renyi) entanglement entropy by subtracting S^ G from S^ EX , 

A Sf = sf EX - sf G . (7) 

We will study the time evolution of Ain following sections. 
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3 Replica Trick 

3.1 Locally Excited States 

In this section, we explain the replica method for locally excited states. The states we 
are considering are given as follows: 

|T) = Me~ iHt e~ eH 0{-l, x) |0). (8) 

Here is a local operator in Shrodinger picture and we introduce the regularization 

factor e~ eH because the norm of the state |T) is divergent without e~ eH and T) does not 
belong to the Hilbert space. This corresponds to smearing the point like excitation. 

The density matrix p is given by 

p = J\f 2 e~ iHt e~ tH 0(—l,x.) |0) (0| 0(—l,x.)^e~ eH e iHt 

= J\f 2 0(r e , —l, x) |0) (0| 0\r h —l, x) (9) 

Here 0(t, —x) is a local operator in Heisenberg picture and in the second line we intro¬ 
duce complex times r e = —e — it and 77 = e — it. In the calculation of (Renyi) entanglement 
entropy, we first treat complex times r e and 77 as if they are real parameters and finally 
we analytically continue to complex values. 

3.2 Replica method 

Now we explain the replica method for locally excited states. In the path integral formal¬ 
ism, we can express the wave functional for locally excited states is given by 

l r<l>(T=Oy)=<p(x 1 ) 

»Wi)) = (zf x ) 3 <0I«> = / 0(T e , -i,x) e - s w. (10) 

J (j)(r=—oo,x l ) 

In the same manner, the we can express the bra vector as follows: 

l r p{t=oo y) 

**(<Kx)) = (Zf x )~ <T|0) = / Vcf) 0 \n,-l,*)e- s W. (11) 

J <j>{T=oy)=<j>(x i ) 

Then, the density matrix for total system is given by the path integral on the space which 
has boundary at r = +0 and r = —0: 

\ptotU' = (0I*> W'> 

= {Zf x Y l [ 2?^t(7i,-I,x)0(r e ,-i,x)e-^, 

( 12 ) 

where Zf x appears in order to keep trp = 1 and it is given by 

p(j)(r=oo,x' 1 ) 

Zf x = / V<f)0\T h -l,x)0(T e ,-l,x)e- s W. ( 13 ) 

J (j)(r=—oo,x' 1 ) 
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Figure. 2: n-sheeted manifold with operator insertion 


Partial trace corresponds to sawing the region which was traced out, so the reduced 
density matrix is given by 

p(j)[r=oo,x' 1 ) 

[Pa X ]h' = {Z^Y 1 / V(t)0 ] {r u -Z,x)£>(r e , -Z,x)e _SM 

X JJ 5(0(+O, x l ) - (j)'(x 1 )) • 5{(f>{-0, x l ) - 0(x‘)). 

x£A 

From this we can see that the only difference between the reduced density matrix for 
ground states and that for locally excited states is the insertion of local operators 0(r e ) 
and O(ri). We need to insert two local operators in each sheet, so finally we need to insert 
2 n local operators in the n-sheeted manifold £ n which is constructed of n flat spaces and 
it has a conical singularity on the entangling surface as in Figure.2. 

Then, the tr(p^ x ) n is given by the partition function with the insertion of 2 n local 
operators: 


t i(p B a X T = (Zf x )~ n [v<pO'(r h 0?)O(r <s ,0?)---O'(r l ,el)O(r e ,el)e- s M. (15) 

where we introduce the polar coordinate (r, 6) on (r, xi) plane , the region of 6 is given 
by 0 < 9 < 2-1171 and 9^ t — 6\ l + 2n(k — 1), see Figure. 2. 03 is almost the correlation 
function and the only difference is that right hand side is divided by (Zf x ) n , not by 
Z n where Z n is the partition function on n-sheeted manifold £„ without any operator 
insertion. If we consider the difference between Renyi entanglement entropy for excited 
states and that for the ground state, we can find that it is expressed by the correlation 
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function on n-sheeted manifold: 



In this way, we can express the difference of (Renyi) entanglement entropy using the 
correlation function on n-sheeted manifold. 

4 Propagator 

In conformal field theories, (Renyi) entanglement entropy for a ground state is invariant 
under the conformal transformation. In order to preserve its conformal symmetry in the 
replica trick, the action on £ n is given by 

s= j (17) 

J Yj n 

where Y l (i = 0,1, 2, 3) obeys the Clifford algebra ({P, P} = 26^1) and 1 is the identity. 

In this case, £ n is given by the flat space except for the origin. We introduce a polar 
coordinate and this geometry is described by 


ds 2 = dr 2 + r 2 d9 2 + dx 2 + dx\. 
The two point function of ^ and is defined by 


(18) 


- <7'</>o(r, 0, x)i/> t (r 2 ,# 2 ,x 2 )) = S 0l (r,r 2 , 0,0 2 , x, x 2 ). 


(19) 


S a b obeys the equation of motion as follows, 



5(r — r2) 5(9 — 9 2 ) 5 (x — x 2 ) 


S a b can be rewritten by 



( 21 ) 


g is defined by 

g(x, x') = 1 ReH(r, r 2 , 9 , 9 2 , x, x 2 ) + Y {) Y l ImH[;r, r 2 , 9 , 9 2 , x, x 2 ). (22) 

where H(x,x') is given by 



( 23 ) 
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G(r, r 2 ) obeys the following equation of motion, 

d 2 1 d 1 d 2 d 2 9 2 \ 

+ —— + oxo + oxo ) G(r, r 2 , 9, 9 2 , x, x 2 ) = 


dr 2 r Q r r 2 dO 2 ( dx 2 ) 2 ( dx 3 ) 2 


8 (r — r 2 ) 8(9 - 9 2 ) 8 (x - x 2 ) 


Next we consider the boundary condition for It is given by, 

9 + 2ri7r, x) = —^(r, 9 , x). 

Then the boundary condition for two point function of if) is given by 

S ab (r, 9 + 2ml, x) = -S ab (r, 9 , x). 

The boundary condition for the green function G is given by 

G(r , 9 + 2mr, x) = e^~ n+1)n G(r, 9, x). 


( 24 ) 


(25) 


(26) 

(27) 


4.1 Computation of Proapgator 

Let’s compute the propagator G(r,r' ,9,9'). It can be expanded by the eigenfunctions 
V(r, 9, x) which are defined by 


CV(r, 9,x) 

where C = (jg, + if + 4^ + 

After expanding it, it is given by 


+ 


(k 2 + k 2 )I/(r, 9, x), 

a 2 A 

(ax3) 2 J ■ 


(28) 


G(x, x') = 


2mr (2n) 2 


E dk 


"O o T| - 

dkk 


re + l + 2 i 


(At) Ji 


-n+l + 2i 
2n 


( kr ') 


l= — OO ' 


'0 


fc 2 + k 2 


(29) 


x e 


x ,) 


where J„(x) is the Bessel function of the first kind. We can rewrite this green function as 

in P21 CHI IISU2ZI ESI E5I - it is given by 


G(x, x') 


87r 2 nrr' 


E 


I -n+l + 2Z 
2n 


|*o+i( 


-re+l+2Z 
2 n 


){e-e’) 


sinh t 0 


(30) 


where to is defined by 

r 2 + r' 2 + lx — x '| 2 

cosht o =- ~2rr' -‘ (31) 

1 By using the map in m, entanglement entropy for the spherical subsystem is map to the thermal 

entropy in the hyperbolic space. This boundary condition in (125(1 agrees with the one which is imposed 

on fermionic field along S 1 direction. 
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If n is odd, the green function is given by 


G(x,x') = 


sinh U 

n 


87 r 2 nrr' s inh t 0 (cosh (^) — cos 
In this case, the boundary condition (1271) is given by 

G(r, 6 + 2 mr, x) = G(r , 0, x). 

The green function given by (1321) obeys this boundary condition. 
If n is even, the green function is 


G(x,x') = 


cos (^f) Sinh 


to 
2 n 


47 r 2 nrr' sinh t {) (cosh ((^) — cos (^yy-)) 
In this case, boundary condition given by fl27|i is given by 

G(r, 6 * + 2 ri 7 r, x) = — G(r , 0, x). 


(32) 


(33) 


(34) 


(35) 


4.2 Analytic Continuation to Real Time 

Up to here, we consider propagators in the Euclidean space. Two local operators are 
located as in Figure.3. We would like to study the time evolution of (Renyi) entanglement 
entropy. Therefore we perform an analytic continuation to real time as follows, 

ti = e- it , 

T e = -e- it , 

where in lorentzian spacetime e is a cutoff parameter which regulates the divergence when 
a local operator contacts with another. 

After performing it, parameters in Euclidean space are related to those in lorentzian 
spacetime as follows, 


(36) 


r 2 — / 2 — £ 2 + e 2 — 2 iet, 
rl = l 2 — t 2 + e 2 + 2 iet, 


cos (6 — 62 ) 
sin {6 — 9 2 ) 


l 2 — e 2 — t 2 

V 2 
2 et 

~V 2 ’ 


V 2 = yj (l 2 — t 2 + e 2 ) 2 + 4e 2 t 2 . 
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(37) 










Figure. 3: The location of operators in Euclidean space. 


4.3 Dominant Propagators 

After performing the analytic continuation in (l36lh we take the limit e —> 0. A few prop¬ 
agators dominantly contribute to n point functions. We call them dominant propagators. 
They are 0(e~ 3 ). If t < l, dominant propagators on are given by 


S a b(r,r 2 ,d - 0 2 ) = 


[«ri + fry, 


167r 2 \// 2 — t 2 e 3 

s “‘ (r ' = lftrV p - w [ltC + (ry ' 

On the other hand, If t > l, they are given by 

1 


(38) 


Sab(r,r 2 ,0- 0 2 ) = 


[«ri + try, 


1 6 tt 2 \/ t 2 — / 2 e 3 

S ab (r,r 2 ,g 2 - 6) = „ A,. + <ry , 


(39) 


167 r 2 \/t 2 — l 2 e 3 

If t <1, dominant propagators on E ra>1 are red arrows in Figure.4. They are given by 

1 =- 3 +fry, 

(40) 


S a b{r,r 2 ,0 - 0 2 ) = 


lQ^VW^¥e 3 1 a ° 

Sabir 2 , r, 0 2 - 0) = ^ „ Jr — [^r° 6 + IT 1 ^] . 


167r 2 i// 2 — t 2 e 3 
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If t > l, they are red arrows and blue arrows in Figure.4. They are given by 


S a b(r,r 2 ,d - 0 2 ) = 


S a b{r 2 ,r,d 2 -9) = 


it + iy 


S a b(r,r 2 ,9 - d 2 + 2n) = 
Sab{r 2 ,r 1 9 2 - 9 - 2n) = 


64ir 2 t\/t 2 — Z 2 e 3 

-(t + iy 

647 i 2 t\/t 2 — l 2 e 3 ab 

it-If 


^r° + —_-r 1 

ab ' , 1 ab 


Q 2 t—l, 

*r°„ + r 1 


ab 


647 \ 2 t\Jt 2 — l 2 e 3 
-it-If 


iT 


2t + / i 


ab 


ab 


r 

t 

° ^ + fp 1 


?T U — 


Sabir , r 2 , 9 -9 2 -2(n- l)vr) = 
S a bir 2 ,r, 9 2 -9 + 2(n- 1)tt) = 


64t rHVt^Pe 3 C ab 

-it-If 


iT 


2f + Z i 


647T 2 tV^ r ^e 3 r ab 
it ~ If ; 

64ir 2 ty/t 2 — l 2 e 3 


t 


r 

o + / T ^ 1 


ab 


?T U — 

11 ab 


ab 


( 41 ) 



Figure. 4: The schematic description of dominant propagators. S a bir, r 2 ,9 — 9 2 ) and 
S a bi r 2, r ,92 ~ 9) correspond to the right arrow and the left arrow respectively. 


4.3.1 Propagators in the Cartesian Coordinate 

Given locally excited states are defined by being acted by local operators on the ground 
state in Cartesian coordinate. Therefore we explain the relation between propagators in 
m and them in the Cartesian coordinate. In the Cartesian coordinate, the action in 
free massless fermionic theory is given by 

S = [ , (42) 


where y^ = {y°, 7 1 , y 2 , y 3 }. 
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When we introduce a polar coordinate x 1 = rcos6,x° = r sin 9 and redefine the 

7 1 7 0 

fermionic field by tj>'(r,9 ) = e 2 u ijj(r } 9), it is rewritten by 


S — dr d6r / d 2 xt/> 


r) r° 

r^ + r 1 — + — + r x a 

r 2 r 




where T° = 7 1 , T 1 = 7 0 , T x = 7 X 

tf(r,0) = cos (^Ji/j(r,d) + sin Q^7°7VM)- 


After perfprming this map, the boundary condition for if ; jj is mapped to 

ip'(r, 9 + 2rar) = (— l) n+1 ip'{r, 6). 


(43) 


(44) 


(45) 


We define propagators by 

- <T^(r,0)^(r 2 ,0 2 )) = S' fe (r,r 2 ,0 - 0 2 ), 

- ?(r 2 ,0 2 )) = S' ac (r,r 2l 0 — 0 2 ) 7 £, = V^r, r 2 . 0 - 0 2 ). (46) 

After performing this map, we also perform the analytic continuation in (j36j) and take 
the limit e —>■ 0. After that, only a few propagators can contribute to n point function 
dominantly. They are 0(e~ 3 ). In any time, dominant propagators on £1 are given by 


S’ ab {r,r 2 ,9-9 2 ) = T ^ 1 t ab , 
S' ab (r 2 ,r,9 2 -9) = ^^ 1 t ab , 
V ab (r,r 2 ,9-0 2 ) = j^l ab , 
V ab (r 2 ,r,9 2 -9) = j^ l ab , 


(47) 


where 1 is identity. 7 0 is changed to i'y* ((7*) 2 = — 1). 

On the other hand, if t < l those on £ ra> i are the same as in (j4T|) . In the region t > l, 

2 When n = 1, i/j 1 is regular at the origin. 
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they are given by 


s'jr^e - e 2 ) = 

S' ab (r 2 ,r,0 2 - 0) = 


— i(t + /) 
647 r 2 te 3 
i(t + /) 
647r 2 te 3 


t — l i , 

—j—lab + 2 7 a& 

1 l llb + 2 7 ab 


t 


S ab( r , r 2,0- 02 + 2tt) = 

S' ab (r 2 ,r,6 2 -6 -2 ti) = 


i(t — l ) 
647r 2 te 3 
i(t — l ) 
647r 2 te 3 


t + l i 0 t 

—j~lab ~ 2 lab 

t + / i „ f 
lab ~ 2 lab 


t 


S'ab( r , r 2 , 0 - 02 - 2 (n » 1 )tt) = (-1) 

S' b (r 2 , r,8 2 -9 + 2 (n - 1 )tt) = (-1)’ 


647r 2 te 3 
64n 2 te 3 


t + l i n t 

— lab ~ 2 lab 
t l I , 

' 7 aft - 2 7 at> 


and 


V a b{r,r 2 ,0 - 0 2 ) = 
V a b(r 2l r,0 2 -8) = 


(t + 0 

647r 2 te 3 

(t + Z) 


t-Z 


647r 2 te 3 


t 

t-l 


bVL+2i- 


t 


(7V) a6 + 2 U 


V a b{r,r 2 ,0~ 0 2 + 2tt) = 
V a b(r 2 ,r,0 2 - 8 - 27t) = 


(*-0 

647r 2 te 3 

(*-0 

647r 2 te 3 


t +1 


t 

t +1 


bVL - 2i- 


OV), 


t v ' ' /ab 


21 


ab 


V ab (r, r 2 , 0 - 0 2 - 2(n - l)7r) = (-l) r 
K 6 (r 2 , r, 0 2 - 0 + 2(n - 1 )tt) = (-l) r 


.1 ft - 0 

647r 2 te 3 
, - 0 


t +1 


(7V) aft- 21 ** 


647r 2 te 3 


t + / 


(7V) 


ab 


21 


ab 


Let’s study the time evolution of AS^ in the next section. 


( 48 ) 


(49) 


5 AS)^ for Various Local Operators 

In this section, we study the time evolution of AS^ l> for various operators by the replica 
trick. Especially, we focus on their behavior in the late time region (i > 1). 

5.1 for i/j r a 

The locally exited state is given by 

1^) = VV>'(-t,-Z,x) |0). (50) 

The time evolution of Ais as follows. 
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5.1.1 The Excess of Renyi Entanglement Entropies 

AS^ for it is given by 


(n) _ 1 


Asr = 


1 — n 


log 


W ( 9 + 2 ( n - l )^Wa( 6 2 + 2 (n - 1 


( 51 ) 


When we take the limit e —* 0, Avanishes in the early time region (t < l). When t 
is greater than l , two diagram in Figure.5 dominantly contribute to AS^ ■ In this region 
the denominator is given by 


(^a(r,9)^ a (r2,d 2 ))^ i 


167r 2 e 3 


(52) 


Therefore Ais given by 


A^ n) = 


Ai — 
A 2 = 


1 — n 

t + r n 
4 1 

t — 1 


l°g [A\ + A 2 ], 
t-l 


t 


(7 7 )aa + 2 , 


(53) 


4t 


t + l 


1\ 


(7 7 


where ( 7 * 7 * ) a a is real because 7 t y 1 is a hermitian matrix^. If we take the late time limit 
(t > 1), Ais given by 


log [Ai + A 2 ], 


A S{ n) = -- 

A 1 — n 

'(7*7%+ 2 


A\ — 
A 2 = 


4 

-(7V)aa + 2 


(54) 


where as in Figure.5, A\ and A 2 are respectively given by 

( v aa (r 2 ,r,e 2 -e)) n 


Ai — 


Ao — 


V»a M)^'(r 2 ,0 2 ) 


Si 


(_l) n+1 (I4 a (r 2 , r, e 2 - e - 2(/r)) n F fla (r 2 , r, 0 2 - 6 + 2(w - 1 )tt) 

^'a{ r i e W a {r 2 ,0 2 ) 


( 55 ) 


3 -l < (7*7^aa < I- 
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If we take the Von Neumann limit n —> 1, A Sa is given by 

A S A = ^((tV )«» - 2 ) log(2 - (7 t 7 1 )aa) - ^((bV)™ + 2) log^yV ) aa + 2) + log(4). (56) 

If (7*7 1 )aa vanishes, A Sa is given by log 2 which is entanglement entropy for the EPR 
state. The lower value of A Sa in (155]) is given by log4/34 (( 7 * 7 ’) a « = 1 or —1). 


5.1.2 Reduced Density Matrix 


If we identify A\ and A 2 as the diagonal components of p n A respectively, it is expected 
that the density matrix for this state is given by 


PA 


1 (hV)aa + 2 0 

4 y 0 -(7VU + 2 


(57) 


where tr aPa and each diagonal components of p A are positive. 

It is expected that the density matrix for T) = A |0) is given by 


Pa 


1 (h W)aa + 2 0 

4 \ 0 -(7VU + 2 


(58) 



Ai A 2 

Figure. 5: The schematic description of Ai and A 2 . 


5.2 Afor '0'^' 

A locally excited state is given by 

|T) = Affi^'i-t, -Z,x) |0), (59) 

which is invariant for SL( 2, C) transformation. 

Let’s study the time evolution of A S A \ 
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5.2.1 The Excess of (Renyi) Entanglement Entropy 

Afor (j59|) is given by 

{^{9 + 27r)-0-0(^2 + ^)i ) '^{9)^(62]) 


A^ 2) = -log 


{^{6)^(6^ 


El 


/ o\ 

In the early time region (t < /), A S\ vanishes when we take the limit e —» 0. 

t > 

A S A ■* = — log [ A\ + A 2 + A 3 ], 


other hand, in the region t > /, A S\ is given by 


Ai = 


Ao — 


c(t + /) 4 

\( 2t ~ i y ii 

2 4 c t 2 (t 2 - l 2 ) 

[\ t) : jj 

c(t - l ) 4 

\f 2t + l \ 2 Jl 

2 4 c t 2 (t 2 - l 2 ) 

1 

h 

1 _ 

1 _ 


2 


A 3 = 2 c 


t 2 -i 2 ^ 2 


2 H 2 c 


5t 2 — l 2 ^ 2 
t 2 
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where c = 4. 


co'i 

If we take the limit t —>• 00 , AS a is given by 

A 4 2) = log 


77 


If t < /, for arbitrary n, AS^ vanishes. In the region t > l, A S^ is given by 


a s { r ] = 


1 — n 


l°g [Ai + A 2 + A 3 ], 


where A 1 , A 2 , A 3 are given by 



A 3 — 2 c 


2 4 ct 2 (t 2 + l 2 ) 
't 2 -l 2 \ n 


2 H 2 c 



t 2 


If we take the late time limit t —> 00 , AS 1 ^ is given by 

n k<n 


a s ( r } = 


1 — n 


log 


12 

64 


+ ' 


64 


22 - c ^ k5 


n—k 


ke 2Z 


(60) 
On the 

(61) 

(62) 

(63) 

(64) 

(65) 
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5.2.2 Reduced Density Matrix 

By using results up to here, we are able to guess the reduce density matrix for |T) = 
as follows. In the limit e —>■ 0, only four diagrams in Figure.6 can contribute to A S^\ 
We define Ax, A 2 , A by 


A 1 = 

a 2 = 

A = 


12 

64’ 

12 

64’ 

if 

64 V 


5 



( 66 ) 


In the replica trick, (Ai) n , (A 2 ) n and 8 • tr(A n ) respectively correspond to the red-lined 
diagram, blue lined diagram and green lined diagrams in Figure.6. Therefore it is expected 
that the reduced density matrix for this state is given by the 10 x 10 matrix, 


PA = 777 


1 

64 


[12 0 0 0 0 0 \ 

0 A 0 0 0 0 

0 A 0 0 0 

0 0 A 0 0 


0 
0 
0 
\0 


0 0 0 A 0 
0 0 0 0 12 / 


where A is given by the 2x2 matrix, 


A = 


5 4 
4 5 


If we diagonalize the matrix A, it is given by 

"9 0 


A = 


0 1 


For any n, (Renyi) entanglement entropy is given by 


A s{ n) = 


1 — n 


log 


2 - 12” + 4- 9” +4 
2 6 " 


(67) 


( 68 ) 


(69) 


(70) 


If we take the von Neumann entropy limit n —> 1, entanglement entropy is given by 

3, /128 \ 

ASU= 4 log b“} 

Min entropy is given by 


( 71 ) 


y(00) 


[! 6 \ 








(Renyi) entanglement entropy for this state monotonically decreases when the replica 
number n increases. Therefore we can consider AS 1 ^ 00 ^ as the lower bound of (Renyi) 
entanglement entropies. 

If we take the large N limit in the U ( N ) or SU ( N ) free massless fermionic held theory, 
the number of diagrams which can dominantly contribute to Adecreases. Because the 
number of trace in the green-lined diagram is less than that in the others, the blue-lined 
and red-lined diagram dominantly contribute to Ain the large N limit. Therefore it 
is given by 


A g(n>2) 


2n - 1 
n — 1 


log 2 + 


n 


n 


log 



(73) 


After that, if we take n —> oo limit, Min entropy is given by 


A S<T> = log 



(74) 


This result describes that Min entropy is same as (1721) even if we take the large N limit. 
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Figure. 6: The schematic explanation of the correspondence between diagrams and com¬ 
ponents of a reduce density matrix. 


5.3 Afor ijj'hjj' 

A locally excited state is given by 

|T) = A-/,x) |0), 
which is variant under the SL( 2, C) transformation. 

5.3.1 The Excess of (Renyi) Entanglement Entropy 

Afor the state in (175]) is given by 

(2) r {^^{e + 27r)^ /t '0 / (^2 + 27r)^'ty / (0)^'V(0 2 )} E 

AS \ — — log -^- 

[ (v*' t </>'(6 | )i/>' , </>'(02))!, 


(75) 


(76) 
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( 2^ 

If we take the limit e —> 0, A S y A vanishes in the early time region (t < l). 
region t > l, Anontrivially grows as follows, 

A S A ^ = — log [Ai + A 2 + A 3 + A 4 ], 






/q\ 

If we take the late time limit (i > 1), A S A is given by 


A4 2) = -log 


20 2 20 2 4 • 3 2 4-3 
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+ -7T + 


+ 


46 46 46 46 

If £ < l, for arbitrary n, A S A ^ vanishes. In the region t > l, A S A ^ is given by 


&sf = - - 

1 — n 


l°g [Ai + A 2 + A 3 ], 


where A 4 , A 2 , A 3 are given by 



+ 4 


+ 4 




If we take the late time limit t —> 00, Ais given by 


Al\ k (3 ^ 2 + / 2 ^ n - fc ' 



f 2 


A S ( A n) = 


1 — n 


log 


/20\ / 3 N 

2 (m) + (m. 


In the 


(77) 


(78) 


(79) 


(80) 


( 81 ) 
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5.4 Reduced Density Matrix 

We define 5 1; 5 2 , 5 by 

( 82 ) 

0 
3 

In the replica trick, (5i) n , ( B 2 ) n and 8 • tr(B n ) respectively correspond to the red 
dashed lined diagram, blue dashed lined diagram and green dashed lined diagrams in 
Figure.7. Therefore it is expected that the reduced density matrix for this state is given 
by the 10 x 10 matrix, 



fll = “ 
64’ 


B, = 


20 

64’ 


1 

B = — . n 
64 VO 


Pa = 


1 

26 


/20 0 0 0 0 0 \ 

0 5 0 0 0 0 

0 0 5 0 0 0 

0 0 0 5 0 0 

0 0 0 0 5 0 

Vo 0 0 0 0 20/ 


(83) 


where 5 is given by 2 x 2 matrix, 


5 = 


3 0 
0 3 


(84) 


by 


If we take Von Neumann entropy limit n —> 1 in (18TT) entanglement entropy is given 


AS A = log 16 + ^ log 2 — | log 3 — | log 5. 

4 o o 


(85) 


Although the local operator ifi'ty' is not invariant under SL( 2, C) transformation, AS 1 / 1 
is invariant under this transformation. 

If we take the limit n —> oo, it is given by 


AS'” 1 = log 


( 86 ) 


If we take the large N limit in the U ( N) or SU ( N ) free massless fermionic held theory, 
the number of diagrams which can contribute to A S^ decreases. Because the number of 
trace in green dashed lined diagram is less than that in the others, blue dashed lined and 
red dashed lined diagram dominantly contribute to A S^ 1 in the large N limit. Therefore 
it is given by 

‘ (n) “ "j ~ log 2 + —^ — log ( -4 
1 — n 1 — n \16 


A ST = 


(87) 
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If we take n —> oo limit, Ais given by 

AS!”’ = log (j ) ■ (88) 

Even if we take the large N limit, the lower bound of Aagrees with that in (186|) 
similarly to the result in the previous subsection. 
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Figure. 7: The schematic explanation of the correspondence between diagrams and com¬ 
ponents of a reduce density matrix. 


5.5 Spin dependence 

Here we discuss the spin dependence of the reduced density matrices in (1571) and ()58lh 
Their components have (7*7' ) 00 . Therefore they depend on the direction of the spin. 
These spin dependences can be understood as follows. To explain why matrix components 
depend on the spin, we diagonalize 7*7* and it is given by 


/I 0 0 0 \ 

0-100 
0 0 10 
\0 0 0 - 1 / 


(89) 


We can derive (po + 7*7 l Pi) *&(p) — 0 from the equation of motion for $(p). The modes 
propagating along x 1 direction (pi 7^ 0,p2 = 0,p3 = 0) obey 


(Po + mV) $(p) 


/Po T Pi 0 0 

0 Po Pi 0 
0 0 Po + Pi 

\ 0 0 0 


0 \ 
0 
0 

P 0 -P 1 J 


$2 (?) 
<h+(p) 

\$4 ip)/ 


0. (90) 


The equation in (1901) describes that that (anti-)particles can propagate along only the 
left or right direction as follow. An anti-particle <f> + (p) (7 t 7 1 $ + (p) = <3> + (p)) can not 
propagate in the right direction parallel to x 1 axis (ad > 0). On the other hand, the 
anti-particle <f> _ (p) (7*7 1 $ _ (p) = —$ _ (p)) is not able to propagate in the left direction 
(ad < 0). Therefore if we define a locally excited state by acting a component of il> or -0 on 
the ground state, their reduced density matrices depend on the representation of gamma 
matrices as in (157)) and (158|) . For example, we choose the basis which diagonalizes 7 t y 1 as 
in (jgtjj) and ipi acts on the ground state. This local operator creates an anti-particle and 
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it propagates with time. At the late time (t > 1), it is necessarily included in the region 
A or B. As explained above, it is not able to propagate in the right direction parallel to 
x l axis. Therefore it is included in A (B) with the probability which is bigger than | (less 
than I). The reduced density matrix in (1F7|) is given by 

* = i (o !) • < 91 > 

The result in f|9T]) describes that the antiparticle which created by 0i are included in A 
(B) with probability | (| ffl If we choose a basis in which (7*7* )n is given by 0, A 
is given by log 2. 


6 Quasi-Particle Interpretation 

In [13 US, US], we have interpreted the late time value of Ain free massless scalar 
held theories in terms of entanglement between quasi-particles. 

In free fermionic field theory, it is expected that the late time values of Acan be 
interpreted in terms of entanglement between quasi-particles. Therefore we decompose 
local operators into left moving modes and right moving modes, 


0a = 0a f + 0^ + 0a + 0a , 

0l = 0a L + 0f + 0a t + 0f, 

0a = 0a f + 0f + i (0V) a + i (0V) O , 

where 0^1 is dehned by 

</7 = • (TV),,, 

^ Rt = i(^» R V) 0 . 


(92) 


(93) 


The left moving mode and right moving mode of the anti-particle and particle are 
dehned as follows. We decompose 0,0t into the momentum modes. Their right moving 
modes (left moving modes) are dehned by the sum of the momentum modes whose ^ 
is positive (^ is negative). As we explained in the previous section, the number of 
momentum modes which the right moving modes (left moving modes) have depends on the 
choice of spin’s direction. Therefore instead of ordinary anti-commutation relationships 
= SRK'dab), we impose the following exotic anti-commutation relationship on 
the particle and anti-particle, 

4 If we act ipi + 02 on the ground state, the late time value of A00 is given by log 2 as we expected 
5 If you choose the bases which diagonalize in 2d massless free fermionic field theory, a component 
of 0 (0) is able to propagate in the only left or right direction parallel to x 1 . Therefore A00 vanish for 
the locally excited states generated by acting a component of 0 (0) on the ground state. 
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{^f} = (tab- \ (tS' 1 ),^ , 
{4>aA L b ] }= ^a6 + ^(7V) a& ) , 
{^a\^b} = (^“^(7V) a6 ) , 
W f »^}= ^a6 + ^(7V) a6 ) , 
{$*$?*} = (^ab- \ (7 1 f) ab ) > 
^a6 + ^(7 1 7 t ) a6 ) • 

The vacuum |0) = |0) L <E) |0) R is defined by 

< 10 } = V-f |0) = fit 10) = o, 


( 94 ) 


(95) 


where K = L, R. The number of momentum modes which the left and right moving modes 
have depends on the choice of 7 t y 1 (the choice of the spin’s direction). Therefore, it is 
reasonable that the anti-commutation for them is given by that in (l94|M We claim that 
the exotic anti-commutation relationship in (|94|) can be applied to Ain 4 dimensional 
free massless fermionic field theory. However in d( 7+ 4) dimensional spacetime, that in (1941) 
should be deformed. For example, in 2 dimensional spacetime the reduced density matrix 
for the state excited by a local operator constructed of chiral and anti-chiral operators 
should become 


Pa 


1(1+ (7* 7 1 )aa 0 \ 

2 \ o 1 - (tVW ' 


(96) 


Therefore the anti-commutation relation in ()94l) should change to another one. 

This way, we can construct reduced density matrices which agree with those which 
are obtained by the replica trick as we will see later. In following subsections we will 
compute them under this decomposition in various examples in 4 dimensional free massless 
fermionic field theory. 


6.1 Various Examples 

Here we compute reduced density matrices for various locally excited states under the 
quasi-particle interpretation. 

6 However we find this exotic anti-commutation relationship heuristically. 
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6.1.1 p A for 0 a 

A locally excited states is given by 


|*>=^a|0>, 


(97) 


which is variant under the SL( 2, C) transformation. Under the decomposition in (l92]h 
this state is given by 


I*) = W l°>L ® I°)r + |0>L ® |0> J , 


(98) 


after the normalization constant A f is determined. 

Normalized orthogonal left and right moving states are given by 



(99) 


where ( M L ) 2 and (A f R ) 2 are given by 



( 100 ) 


A reduced density matrix pa is defined by tr^p. It is given by 


PA = \[( 2 + (7V) J |0> fi (0| H + (2 - ( 7 V) J \tf) R WIJ (101) 


which agrees with the one obtained by the replica trick. When (7*7 1 ) aa vanishes, A 


are given by log 2. In the n —> 1 limit, it is given by entanglement entropy for maximally 
entangled state (EPR state). If we choose gamma matrices as in (j89lh the anti-particles 
created by 'ipi is not able to propagate in the right direction parallel to x 1 axis. Therefore 
the probability with which it is finally included in B is larger than the one with which it 
is included in A at the late time. The reduced density matrix in (1101(1 agrees with our 
intuitive expectation. 

6.1.2 pa for '00 

A locally excited states is given by 


|T) = A70-0 |0), 
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( 102 ) 








which is invariant under the SL(2, C) transformation. 

By using an unitary transformation, y^ 1 can be diagonalized. After diagonalize it, 
the reduced density matrix for this state can be computed easily. Here it is given by that 
in (159]) 

Under the decomposition in (I92[) . the state in (I102[) is given by 

in = \ £ iw+evf+y 8 w+wi io) L ® io> K , (103) 

a 

after the normalization constant A f is determined. 

Normalized orthogonal states are given by 

|0>L, 

| = loy. 

a 

|^V% = N RR |0>* , (104) 

W) =-^'a L V'a t | 0 ) L , 

W ) = 0) L , 

m=KV\0)n, 


where normalization factors are given by 


(a r LL ) 2 = 

: (a r RR ) 2 = 

1 

“ 3’ 







(A?) 2 = 

( A *) 2 = 

«- 


= (Aff) 2 

= m 2 = 

- (m 

1 = m 

2 = 2, 

( a ?) 2 = 

= 

M) 2 ~ 

M) 2 

= W 8 ) 2 

= (A?) 2 = 

- m 2 

' = (AC 4 8 ) 

2 _ 2 
“ 3' 









(105) 

A reduce 

density matrix is given by 






= tr L p = 









12 . . 

“ 64 |0)iJ 

(°l R + 

U R ) 

64 'R 

W\r 

64 

A 8 }* AH 

12 

H+ 64 

i*v% 



(106) 


where Mi and Mi are given by 

Adi = M.s = M.2 = M-i — 1, 

Ad 2 = Ad 4 — Adi = Ms — 9. 

This reduced density matrix agrees with the one obtained by the replica trick. 


(107) 
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6.1.3 pa for ipty 

A locally excited states is given by 


|T) = A |0), 


( 108 ) 


which is variant under the SL( 2, C) transformation. However Afor (1108j) is invariant 
in the replica trick. Therefore we diagonalize 7 t y 1 as in (l89]h After diagonalizing it, pa 
is given by 


p A = tr R p = H (<t> R ^ R \ R + ^ |0) 


64 


R \ u l R 


h Y W) R Wl r + ^zY W) R 


R I 

* \R ‘ 


( 109 ) 


7 Conclusions and Discussions 

In this paper we have studied the dynamics of quantum entanglement. We considered the 
4 dimensional free massless fermionic field theory. Locally excited states are generated 
by being acted on the ground state by various local operators. We defined the excesses 
of (Renyi) entanglement entropies Aby subtracting (Renyi) entanglement entropies 
for the ground state from those for locally excited states. Aare finite quantities. We 
found that Avanish if t < l. If t > /, Aincrease and finally approach constants 
as in scalar held theory mmm- 

We found that if the locally excited state is given by Afijj a |0) ,Afip a |0), components of 
its reduced density matrices include 7 t y 1 . Therefore they depend on the direction of spin. 
A component of ^ or ^ creates an anti-particle or particle. Their propagation depends 
on the direction of spin. Therefore components of reduced density matrices for those 
locally excited states depend on the direction of spin. We also generated locally excited 
states AT^ip |0) and A |0) and investigated the time evolution of Afor them. We 
found that Aare invariant under the lorentz transformation although A A/4-0 |0) is not 
invariant under the lorentz transformation. 

We found that the time evolution and the late time values of Acan be interpreted 
in terms of quasi-particles. However in the fermionic held theory, the number of modes 
which propagate in the left ( x 1 < 0) or the right (x 1 > 0) direction depend on the choice of 
the spin’s direction. Therefore we need to impose exotic anti-commutation relation in (l94j) 
on quasi-particles. In the 4 dimensional free massless fermionic held theory, the results by 
the replica trick can be interpreted as the quantum entanglement between quasi-particles 
if the relationship in (1941) is applied to them. In the d(^ 4) dimensional free fermionic 
held theory, anti-commutation relationship should be changed. 

We would like to list a few of future problems: 

- It is interesting to study charged (Renyi) entanglement entropies for locally excited 
stats and holographic held theory. 
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- It is important to study how we should deform the relationship in (1941) in d{^ 4) 
dimensional free massless fermionic field theory. 

- It is interesting to study the time evolution of AS^ 1 ' 1 in non-relativistic field theory. 

- It would be expected that gamma matrices included in reduced density matrices 
is related to the modular Hamiltonian. It is interesting to clarify the relationship 
between the spin dependence of the reduced matrices and the entanglement Hamil¬ 
tonian. 
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